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Abstract 

We study Smoluchowski-Poisson equation in two space dimensions provided with Dirichlet 
boundary condition for the Poisson part. For this equation several profiles of blowup solution 
have been noticed. Here we show the residual vanishing. 


1 Introduction 

We study parabolic-elliptic system proposed in statistical physics to describe the motion of mean 
field of many self-gravitating Brownian particles [2D]. It is composed of the Smoluchowski part 

ut = A u — V • uVv in 9 x (0, T) (1) 

with null-flux boundary condition 


du 

dv 



du 

- u— = 0 
ov 

on (9fl x (0, T) 

(2) 

and the Poisson part in the form of 





— Av = u, 

o 

II 

G 

CO 

(3) 


where U C R 2 is a bounded domain with smooth boundary dQ and v is the outer unit normal 
vector. Initial condition is given as 


u \ t=0 = uq(x) > 0 in tl, (4) 

where uq = uo(x) is a smooth function. 

System m-m is subject to thermodynamical laws, total mass conservation and free energy 
decreasing, 


d f f fdudv 

— u = /V ■ (Vm — uVv) = —- u— ds = 0 (5) 

dt Jn Jn Jon ov ov 

— F{u) = — / u|V(logu — v )\ 2 < 0 (6) 

dt Jft 
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where ds denotes the surface element and 


F{u) = f u(logu — 1) dx — tt((—A) 1 u,u) (7) 

in 2 

with v = (—A)~ 1 u standing for Q. 

A related model arises in the context of chemotaxis in theoretical biology IS S], where the 
Poisson part is provided with the Neumann boundary condition such as 


a 1 f dv 

V ^ U ~W\Jn U ’ to 


= 0 , 


on 



( 8 ) 


Concerning (jU)-d2]) , dU), and ©, there is a threshold of ||uo||i = A for the blowup of the solution. 
More precisely, if A < 4-7T the solution exists global-in-time mmm- If a local mass greater 
than 47T is concentrated on a boundary point, on the contrary, there arises blowup in finite time 
Mm- Underlying blowup mechanisms were suspected from the study of stationary solutions [2], 
This attempt was followed by mm, using radially symmetric and general stationary solutions, 
respectively. Up to now several properties have been known nanaiE]. See our previous work [21] 
and the references therein. System dU)-([4]) , provided with Dirichlet condition for the Poisson part, 
is taken by [23]. It excludes boundary blowup points. Here we continue the study [2][j on interior 
blowup points. 

Fundamental features of system m-® are the following. First, local-in-time unique existence 
of the classical solution is standard, given smooth initial value uq = uq(x ) > 0. Henceforth, 
T G (0, Too] denotes its maximal existence time. If uq ^ 0, which we always assume, the strong 
maximum principle and the Hopf lemma guarantee u(-,t ) > 0 on for t > 0. Maximal existence 
time T of non-stationary solution u = on the other hand, is estimated from below by Huolloo- 

Hence T < Too implies 

lim ||u(-,t )|| 00 = Too 

t\T 

and the blowup set 


S = {xq G | 3xk —>■ xo, 3tk t T such that u(xk,tk ) —>• Too} (9) 

is non-empty. Since boundary blowup points are excluded [23], if T < Too in (P)-(ll]) we have 

u(x,t)dx —*• ^2 m ( x o)fix 0 (dx) T f(x)dx in A4(fl) = C(Q)' (10) 

xo e<S 

as t t T. Here, the blowup set satisfies 5 C SI with ([5 < Too, and it holds that 0 < / = f(x) € 
L 1 (0)nC(U\5). 

The blowup mechanism at each inner blowup point, described by [24], is more complicated than 
the ones suspected before. Let xo € S and R(t) = (T — t) 1 / 2 . As is shown in [22, 25] it holds that 

lim limsup|||u(-,t)|| L i (B( XOjb m t )) — 7n(x 0 )| = 0. (11) 

6t+°° t\T 

Henceforth, Cj, i = 1,2,--- ,15, denote positive constants. 
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Theorem 1 ([23]). Any tkf T admits a sub-sequence, denoted by the same symbol, and m € Nu{0} , 
provided with the following property. Thus, given 0<£< 1 and R^> 1, there is s > 1 such that, for 
t' k t T defined by R(t ' k ) = sR(t k ) we have 0 < r J k < C\R(t ' k ) and xj, € B{xo,C\R(t' k )), 1 < j < m. 

m 

Then, for B k = B(x k ,bj) and Gk = |^J B k , it holds that 

3 = 1 


B k n B l = 0 , i ^ 3 , k > 1 


lim sup 

k—> oo 






— 87T 


< e, 1 < j < m 


lim limsupi2(t , fc ) 2 ||«(-,t , fc )|| L oo (B( * 0l 6fl ( t' )) \G fc) < Cfe-R 2 . 

oj+oo k^-oo 


( 12 ) 


In this paper we show that if (1121) is extended continuously then it holds that m(x o) € 87 tN. 
Namely, here we assume that any e > 0 admits rrij(t ) € N U {0}, Xj(t) € B(xo,C^R(t)), and 0 < 
rj(t) < C^Rft) for 0 < T — t <C 1 and 1 < j < m(t) € Nu{0}, such that, for Bj(t) = B(xj(t),rj(t)) 

m(t) 

and G(t) = [J Bj(t ) it holds that 
3 =i 


Bi(t) n Bj(t) = 0, i/j 


lim sup max 

tfT 1 <j<m(t) 




< £ 


lim limsupi?(t) 2 ||u(-,t)|| i oc (i3(a . 0!fci?W) \ GW) < C 4 . 
“T+oo t|T 


Then our result arises as follows. 


(13) 


Theorem 2. If U3\) holds in with T < +oo, then there is tkf T such that 

\\<MlHB( X 0 ,m(t k ))\G(t k )) < £ ■ ( 14 ) 

Hence we obtain m(x o) = Sirm, m € N, in GW, and in particular, m{t) = m in 17311 . 

We call (1141) the residual vanishing. Although m(xo) € 87tN follows from (1111) . (1131) . and ED, 
we shall show m{x o) € 87 tN first, and then (1141) . In future we shall discuss the problems to derive 
m from (gZD and also to refine ED to 


lim 

frt+OO 


lim sup 

t\T 


l^("> ^) III/ 1 (B(xo,bR(t))\G(t)) < 


£. 


This paper is composed of three sections. Taking preliminaries in the next section, we prove 
Theorem [2] in the final section. 


3 






2 Prelimiaries 


Weak solution is a fundamental tool in later arguments. This notion was introduced first for the 
pre-scaled Smoluchowski-Poisson equation pj3j. Let G = G(x,x') be the Green’s function to (|3|). 
Then we say that 0 < p = fi(dx,t) E C*([0,T],A4(f2)) is a weak solution to (HI)-© if there is 

0<v = v(;t)eL?(0,T-,£') 


called multiplicate operator satisfying the following properties, where £ is the closure of the linear 
space 

,2/7W , &<P 


£o = {'f’ + P<p | V’ £ C(£l x fi), ip E X}, X = {ip E c (yt) | — 


= 0 } 


an 


in L°°(Q x n) and p v (x,x') = Vip(x) ■ \7 x G(x,x') + V<p(x') ■ V x /G(x, x'): 

1. For ip G X the mapping t E [0,T] (ip, p(dx,t)} is absolutely continuous and there holds 


^(<75 //(da:, f)) = (A<^, p(dx, t)) + *)>£,£' a.e. L 


2. We have 


K-^IccnxH) = M^M) ®p(dx',t) a.e. t. 

Here we confirm that the property ^ > 0 of v E £’ means 

| (/»£,£'| < (s,*') 

for any f,g E £ satisfying |/| < g a.e. in fl x fL 

We note the following properties. First, the total mass conservation of this weak solution 


(15) 

(16) 


p(£l, t) = p(£l, 0), tE[0,T] 

is obvious. Next, this weak solution cannot be a measure-valued solution constructed in PI 0133!. 
In fact, any collision of collapses are not admitted here, and more precisely, we have the following 
property. 

Lemma 2.1 ([15]). If the initial meause po(dx) E M(£l) admits xq E such that 

Po(i x o}) > 87T, lhn-^2 (\ x - x o\ 2 Xb(x 0 ,R), Po{d x )) = 0 

then there is no weak solution to m-rn even local-in-time. 

It is, however, provided with the following property, derived from the fact that £ is separable. 

Lemma 2.2 ([IS]). Let {p k (dx,t)} C (7*([0,T],_M(fl)) he a sequence of weak solutions to (©-(©). 
Let the associated multiplicate operator of p k (dx. t) be E L£°(0, T; £'), and assume 

p k (£l, 0) + sup ||&'fc(-,i)||£' < G 5 , k = 1,2,--- . (17) 

te[o ,T] 
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Then we have a subsequence denoted by the same symbol, /i(dx,t ) G C*([0, T], A4(Q)), and v(-,t) G 
L£°(0, T; £') such that 

^i k (dx,t) —*• [i(dx,t ) in (7*([0, T], .Ad(fl)) 
v k{-,t) ^ v{-,t) in L™(0,T-,£') 

up to a sub-sequence, and this p,(dx,t) is a weak solution to with the multiplicate operator 

v(-,t) satisfying 

/i(O,0) + ||^(-, < C 5 . 


We agree with the following notations. First, if p,(dx, t ) has a density as 
ix{dx,t) = u(x,t)dx, 0 < u(-,t ) € C([0, T], L 1 (n)), 
then the multiplicate operator is always taken as 

v{-,t) = u(x, t)u(x', t) dxdx', 

recalling £ C L°°(£l x fl). Under this agreement, condition (1171) is reduced to 

Ax*(n, 0) < C 6 , fe = 1,2,--- (18) 

if each /j, k (dx,t) takes density in [0,T) such as 

H k {dx,t) = Uk(x,t)dx, 0 < Uk = Uk(-,t ) G C([0, T), L 1 (kl)). 


In fact, since we take 




in this case, inequality (fT8l) means 


u k (x, t)u k {x ', t)dxdx' 


IKM)lli = IKdcO) 0)l|i = 0) — X k < Cq. 


Therefore, (1171) follows with 

\H;t)y = \l<cl 

Using this property, we shall derive a hierarchy of weak solutions in later arguments. 

We can define also the regularity of the weak solution. First, given p = //(•, t) € M(Q), we have 
a unique v = v(-,t ) € W 1,q (£l), 1 < q < 2, such that 


-Au = /x, u| an = 0. 

Let I C (0, T) be an open interval and w C 11 an open set. If the weak solution n(dx, t ) has a 
density u = u(-,t ) € L p (u) in u C fl, 1 < p < oo, for t G /, the above u = u(-,f) is in (u>) 
from the elliptic regularity. By Sobolev’s and Morrey’s imbedding theorems, this regularity implies 
(uVv)(-,t) G L\ oc { oj). Then we assign 

-(<p,IJ,(dx,t)) = (A<p(dx,t))-{-('V(p-'Vv,iJ,(dx,t)), a.e. i € / 
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for any p E Cq(uj). In such a case we say that p(dx, t ) is regular in u x I. In Lemma 12.21 if pk(dx, t ) 
is regular with the density Uk(x,t ) in w x (0, T) satisfying 

SUp \\u k (-,t)\\ L pr u) < C 7 
te[o,T] 

for p > 1, then the generated p(dx, t ) is also regular in cj x (0, T). Conversely, we have the following 
properties by the e regularity mm- First, if the weak solution /j,(dx,t) E C*([0, T], is 

generated by a sequence of classical solutions {uk(x,t)} then its singular part p s (dx,t) is composed 
of a finite sum of delta functions. Furthermore, if p(dx,to), 0 < to < T, is regular in the sense of 
measure in an open set w C O, then it is regular in the above sense. More precisely, p(dx, t ) takes a 
smooth density function / = f(x, t) in u x (to — S, to + 5), where u is an open set satisfying wCii 
and 0 < 5 <C 1. 

The weak solution a(dx,t) € C*(—oo, +oo; .A/f(R 2 )) to 

at = A a — V • aVr * a in R 2 x (—oo, +oo) (19) 

is defined similarly, where T(x) = — log -—- is the fundamental solution to —A, Ad(R 2 ) = Co(R 2 ) / 

27r \x\ 

with 

C 0 (R 2 ) = {/ E C(R 2 \J{oo}) | /(oo) = 0}, 

and R 2 U{°°} denotes one-point compactification of R 2 . Then the following property is shown, see 

Pil- 

Proposition 2.1 (Liouville property). Let 0 < a = a(dx,t ) E C*((—oo,+oo), A1(R 2 )) be a weak 
solution to [19\) with uniformly bounded multiplicate operator. Then we have either a(R 2 , t) = 87r 
or a (R 2 ,t) = 0, exclusively in t & R. 

We can also define the weak solution £(dy,s) E C*(—oo, +oo;7W(R 2 )) to 

Cs = AC-V-CV(r*C + |y| 2 /4) in R 2 x (—oo,+oo), (20) 

which arises as the weak scaling limit of u = u(x,t). Thus, given xo E S C II, we take the backward 
self-similar transformation 

z{y,s) = (T -t)u(x,t), y = (x - x 0 )/(T - t) 1/2 , s = -log(T-t). 


Let tfc f +°° aR d put 


s k = -log(T-f fc ) t+oo- 

Then, passing to a sub-sequence denoted by the same symbol, we have 

z(y,s + Sk)dy ((dy,s) in C*(-oo, +oo; AI(R 2 )), (21) 

where C( dy,s ) is weak solution to (f20l) provided with a uniformly bounded multiplicate operator, 
in m, the important property called parabolic envelope arises as 

C(R 2 ,s) = rn(x 0 ) > 0, (\y\ 2 ,({dy,s)) < C 8 (22) 
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valid to s £ (—oo,+oo) with m(x o) > 0 defined by (ITU]) , see [22 1 [25]. 

Here we take the scaling back of C(dy, s), defined by the transformation 

A(dy', s') = e s ((dy, s ), y' = e~ s/2 y, s' = -e~ s . (23) 

It has an extension as 0 < A = A(dy,s) € C*((—oo,0],.M(R 2 )) with A(dy, 0) = m(xo)5o(dy). It 
becomes also a weak solution to 

A s = AA - V • HVr * A in R 2 x (—oo, 0) (24) 

satisfying 

H(R 2 ,s) = m(x o), — oo < s < 0 (25) 

with a uniformly bounded multiplicate operator. 

Now, given sn f +oo, we take 


An{dy) = A(dy,-sn)/m(x 0 ) 

(26) 

to apply concentration compactness principle [7j (see also p. 39 of [21]). 
following lemma, which implies Theorem [TJ 

Then we obtain the 

Lemma 2.3 (concentration compactness). Passing to a sub-sequence we have m € N U {0} such 
that any e > 0 admits yj E R 2 and bj >0, 1 < j < m, satisfying 

lim \y\ - y\\ = +°o, Vi^j 

a — ^00 

(27) 

lim sup | An(Bj)) - 87t| < e, Vj 

£—KX> 

(28) 

\y{ I < Cg(l + ma JLbj)sV 2 , W > 1, Vj 

7 

(29) 


for Bj = B(yj,bj). Furthermore, there arises one of the following alternatives. 

m 

1. m(x o) > 87 vm + e and An, t 3 > 1, is regular in R 2 \ |^J H(y^, bj). It holds that 

i =1 


2. m(x 0 ) = 8irm and 


lirninf An ( R 2 \ I J Bj ] > m(x 0 ) — 8irm — e 

£—►00 ' ' s - / ' ' 


e= 1 


n^o H^lli' x> (R 2 \U™iSi) “ °’ 


lim sup An 

£—►00 


m \ 

R 2 \ U B(yg,bj)\ 


< £. 


(30) 

(31) 


(32) 
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3 Proof of Theorem [2] 


From the assumption ((dy, s) generated in the previous section satisfies an additional condition. 
Namely, each 0 < e <C 1 admits si 2 > 1 provided with the following properties. First, for s > s± 
there are m(s) £NU {0}, y 3 (s) € R 2 , and bj(s) > 0, 1 < j < m(s), such that \yj(s)\ < Cios 1 / 2 and 
bj(s) < Ciq. Next, ((dy,— logs) is regular in |J (R 2 \ Eg) x { - log 5} for 

S>S1 

m(s) 

Eg = J Bj{S) t Bj(s) = Bis-^y^rs-^bjis)). (33) 

3 = 1 

Finally, it holds that Bi(s) FI Bj(s) = 0, i / j, and 

sup ||C(dy,-loga)||i,oo (R 2 \ £ , s ) < Cn 

S>S 1 

sup |C(Rj(s), - logs) - 87 r| < e. (34) 

s>si, l<j<m(s) 

If m(x o) 0 87 rN, we have always the first alternative in Lemma 12.31 which implies the existence 
of 5 > 0 such that 

inf C(R 2 \ Eg,- log s) > 6. (35) 

S>S 1 

If (1351) is not the case there is S 2 > si such that 

C(R 2 \ E S2 , — log S 2 ) < £■ 

Then it holds that (|14l) with some tk f T. Hence Theorem [2] is reduced to the following lemma. For 
the proof we use (1201) . particularly, the term |y| 2 /4, which attract the density of ((dy, s ) to |y| = 00 . 

Lemma 3.1. The weak scaling limit £( dy,s ) € C*(— 00 ,+ 00 ; A1(R 2 )), generated by H21\) . does not 
satisfy and \35 1) . simultaneously. 

Proof. Similarly to Lemma imi concerning (HI)-©, given §1 f + 00 , we have a sub-sequence denoted 
by the same symbol such that 

((dy, s — §g) C{dy, s) in C*(— 00 , + 00 ; A1(R 2 )). (36) 

This ((dy, s) is a weak solution to d20l) . Furthermore, since {((dy, s — s^)} is tight by (1221) . it holds 
that 

C(R 2 ,s) = m(x 0 ), (\y\ 2 ,((dy,s)) < C 8 . (37) 

We have also 

s~ 1/2 \Vj(s)\ < C w , lim s~ 1/2 bj(s) = 0 

sf+OO 

in (1331) . 

Similarly to the remark after Lemma 12.21 the singular part of ((dy,s ), s € R, denoted by 
( s (dy, s), is composed of a finite sum of delta fucntions. By applying Lemma l2Jl to the scaling back 


A(dy, s ) of ((dy, s ), defined as in (f23l) . the coefficient of each delta function of Cs(dy, s) must be less 
than or equal to 87 t. These properties guarantee that the singular support of ((dy, s ), denoted by S s , 
is composed of a finite number of collisionless accumulating points of {s~ l ^ 2 yj(s) | 1 < j < m(s)} 
defined for — logs = s + §i as i —» oo. We may assume also that S s , s G Q, is composed of their 
converging points by a diagonal argument. 

Therefore, we have m(s) € N U {0} and ylo(s) G R 2 , 1 < j < m(s), such that 

5s = {yL(s) | 1 < j < m(s )}, C({y^o(s)}, s) = 8vr 

for s G R. Furthermore, the density function g = g(y, s ) of the absolutely continuous part of £ (dy, s ) 
is provided with the properties 

0 < g = g(y,s) < Cu, / \y\ 2 g(y,s) dy < C 8 , ||c?(-,s)lli < m(x 0 ) 

J R 2 

Ils(-,s)||i > s, s G Q. (38) 


It holds also that S s C Br for any R > C\q. 

Since £ {dy, s) G C*(—oo, +oo; A4(R 2 )) the set G = U S (R 2 \5 S ) x {s} is open in R 2 x (—oo, +oo). 
Furthermore, the above g = g(y, s ) is smooth in G and it holds that 


g s = A g — V • gVw in G 

2 rn{s) 

w(y, s) = -j- + 4 log 

3 = 1 


\y-yL(s)\ 


+ r * g. 


Let 

v{y, s ) = |y| 2 /4 + (r * g)(y, s ). 

Henceforth, we put 

— 1 r 2n 

/M = 2 ^. J o f( re )d6 

for / = f(y), y G R 2 , where y = re 10 is the polar coordinate. 

We take B r = B(0, r), r > R, and so G R, to set 


(39) 


(40) 

(41) 


m(so) 

B; = B r \ j B(yUso),e). 

3=1 

By Lemma [2.11 again, any collision of the points in S s does not occur as s varies. Therefore, making 
0 < £ < 1, we obtain ft (^B(ylo(so), e) n5 s j < 1 and 


y l oc(s) G £(y2o(s 0 ),e) 


(2/ - y l oo(s)) -u y > 0, y G dBiy^so), e) 


for 1 < j < m(so) and |s — so| 1, where v y denotes the outer unit normal vector. 
Then, (1391) implies 


d_ 

ds 


/ g(y,s)dy< g r (y,s)dcTy- (gv r )(y,s) da y 

JdBS. JdBi. 


(42) 
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for |s — so| <C 1, where da = da y denotes the line element. We integrate (H2j) in t, to convert it to 
an inequality on the difference quotient with the mesh h > 0. Now, making e J. 0 and then h ], 0, 
we obtain 


where 


d + 

ds 


9 dy < 




(gv r ) da 



9r~ ~g~ g {r * g) r da, 


A(s) = limsup -^(H(s + h) — H(s)). 


ds 


hiO 


h 


Then, inequality (j43]) is valid to any r > R and —oo < s < +oo. 
From || 5 (-,s)||i + || 0 (-,s)|| oo < Ci 2 , it follows that 


(43) 


sup||(vr*g)(-,s)||oo < +oo. 


Then, (1431) implies 


d+ f 
' 9 ~ 


B r 

d r2n 


‘ 9r- -9 + Ci3 g da 

dB r 2 


r 2 tt r 2 


r 2-rr 


= —(r / g d6) — / (l + —)gdO + C 13 r g d0, 


dr 


i o 


which means 

r { ^ j. 2 

— / rgdr < —( rg ) - (— - Ci 3 r + 1 )g, r > R, -oo < s < +oo, 
ds J q dr 2 

recalling (141 jl . Here we have 

f r Cl4 

B(r,s) = J rg dr > 6 - 9 ^ ^ , r > R, s€ Q 


(44) 


(45) 


(46) 


by (13811 . The strong maximum principle, on the other hand, implies also B(r, s) > 0 for any (r, s). 
Inequality (|45l) means 

d+B < B rr - a{r)B r , o(r) = ^ - Ci 3 + -. 

2 r 

r 2 

Using A(r) = — — Ci 3 r + logr, we have a = A! and hence 

d+(e~ A B) < (e~ A B r ) r , r > R, —oo < s < +oo 
by B rr — a(r)B r = e A (e~ A B r ) r . Now we take r\ > R such that 

r C~\ a 

a(r)>- + 1, 6- 7 > S/2, r>n. 

4 r z + 1 
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Let 0 < p = <p{r) be a C 1 function on [ri,oo), piecewise C 2 , satisfying <p{r\) = 0, ip(r) > 0 for 
r > ri, and 

roo 

/ e~ A p dr < +oo, lim e~ A (tp + ip r ) = 0. (47) 

Jri rf+oo 

Then we put p{r) = 0 for r € [0,ri]. Since B > 0, B r > 0, B(oo,t) < C 15 , each s € R admits 
Tj f +00 such that B r (rj, s) —> 0, which guarantees 

— / e~ A Bp dr < ( e~ A B r ) r ip dr = — e~ A B r p r dr 

ds J 0 Jo Jo 


< 


roo roo 

/ ( e~ A ip r ) r B dr = ( ip rr — cup r )e~ A B dr. 

Jo Jo 


We impose, furthermore, that the existence of p > 0 such that 

ip rr — aip r < —p<p, r > r\ 

To assure all the above requirements to tp(r), we take 0 < c\ <C 1, for example, and put 


(48) 

(49) 


<p(r) = 


ci?’ + c 2 , r > r 2 

sin j3{r - ri), r\ < r < r 2 = n + ^ 


where 


7T 


r2 = ri + I717’ C2 = 75 (1 -i ) - ri ‘ :i ' ' 3 =^ c - 


Then we see 0 < tp = p(r) € C 1 ^, 00 ), <p(r{) = 0, and (j47]h Making ci | 0, on the other hand, we 
obtain r 2 t + 00 . Therefore, (jUJ) arises for 0 < p <C 1 by 


^Prr — 


0 , r >r 2 

—fJ 2 <p, ri < r < r 2 


and 


r 1 

a(r) = o ~ ^13 + “> <Pr = ci > 0, r > r 2 . 
2 r 


Since (HU) is obtained we have 


dJ_ r 

ds Jo 


e A Bp dr < —p I e A Bp dr , s € R. 


i 


—A ; 


This means 


and hence 


d + 

ds 


roo 

et* I -~ A 


e A Bp dr j < 0, seR 

rOO 

lim / e~ A B(-, s)ip dr = 0. 

st-t-00 J o 


We have, on the other hand, 

roo r 

/ e~ A B(-,s)tp dr > - 


e A <p dr > 0, seQ 


by (1461) . a contradiction. 


□ 
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